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Abstract 



We find the explicit operatorial form of renormalon-type singularities in 
abelian gauge theory. Local operators of dimension six take care of the first 
U.V. renormalon, non local operators are needed for I.R. singularities. In 
the effective lagrangian constructed with these operators non local imag- 
inary parts appearing in the usual perturbative expansion at large orders 
are cancelled. 



I. INTRODUCTION 



It is well known that in field theory the perturbative series diverge and are asymp- 
totic at best so that the question of the reliability of the perturbative expansion is raised. 
In particular the factorial growth of the perturbative coefficients has been related to the 
existence of instanton solutions || and, in this way, a great number of results has been 
found (for a review see 0). The perturbative expansion for renormalizable field theories 
is plagued by further divergences, the so-called renormalons, which show up as singulari- 
ties on the real axis of the Borel transform of the perturbative series Moreover this 
kind of problems has recently gained a great interest fi-[l3[| and many questions, mainly 
concerning the asymptotic behaviour of perturbative expansions in QCD, still wait for an 
answer. 

The standard technique to cope with divergent series is the introduction of the Borel 
transform. If a function F{a) admits an asymptotic expansion and satisfies certain ana- 
lytic properties [13,03], it can be uniquely reconstructed from the Borel transform F(b) 



defined by: 

F(a) = r dbe- b/a F(b) (1.1) 



The presence of singularities (poles or cuts) on the positive real axis of the Borel transform 
forces to give a contour prescription (for example a principal value prescription) in order 
to define the perturbative expansion via the integral ( |1 . 1|) . This results in an ambiguity in 
the sum of the perturbative series given by the integration above and below the positive 
real axis in the Borel plane. 

To be definite, a singularity of the form: 

~ (L2) 

gives rise to an imaginary part in the function F(a) of the form: 

ImF(a) ~ a s e~ bo/a (1.3) 

which in turn is related, via a dispersion relation, to the large order behaviour of the 
coefficients f n of the perturbative series: 

f n ~ const (^-] n s T(n-s) (1.4) 



The renormalon-type singularities occurr for 

b = ±n/p 2 (1.5) 

where n is a positive integer and j3 2 is the first term of the /3-function and the sign plus 
refers to the ultraviolet renormalons while the sign minus to the infrared ones. These 
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divergences are usually the closest ones to the origin and give the leading asymptotics for 
the perturbative series. 

In the following we want to study the structure and a possible classification of the 
singularities due to renormalons, in the framework of QED with a large number of flavours 

g. 

The aim of the first part of the paper is to present a prescription which could allow to 
take sistematically into account the imaginary parts (see eq. ( p..3| )) coming from ultraviolet 
renormalons. The idea stems from a conjecture made for the first time by Parisi for the 
4 theory jnj and explicitly checked in the 1/N limit of this theory [17-19 . 



This conjecture relates the imaginary parts coming from the renormalon at b = n/ftz 
to the insertion at zero momentum of local operators of dimension d = 2n + 4. Renormal- 
ization group arguments fix the form of the imaginary parts from which the asymptotic 
behaviour (|1.4j ) is extracted (with given 6 and s) up to an overall normalization constant. 

We implement this conjecture adding gauge-invariant operators as imaginary coun- 
terterms in the lagrangian. The feasibility of such a description for the renormalon-type 
singularities relies on the fulfilment of two conditions: 

• A unique choice must exist for the operators in order to remove the imaginary parts 
of all Green functions. 

• In the effective theory described by the lagrangian with the counterterms, non local 
singularities in the Borel transformed Green functions must be absent. 



The second half of the paper concerns the problem of infrared renormalons [StKS] 



In the case of QCD a connection has been conjectured between infared renormalons and 
the admissible multilocal gauge-invariant operators which can be seen as the dual of 
the local operators of the theory [p0[| . Then from the absence of local gauge-invariant 
operators of dimension two, the absence of the infrared renormalon at b = —1/(32 has 
been argued. Recently the possibility of the existence of this renormalon has been raised 
[|]|| and this would suggest the presence of a 1/q 2 term in the operator product expansion 
for euclidean correlation functions [^3j and would spoil the connection with the gauge- 
invariant operators. 

While on general grounds this possibility cannot be ruled out, an explicit calculation 
23J shows that this pole doesn't indeed exist within the 1/Nf expansion for QED. More- 



over QCD is believed to have the same behaviour of QED after replacing the right /?2 
coefficient. 



In the spirit of the calculation of p4j we study the structure of the first infrared 
renormalons in QED. In particular we find the form of the multilocal counterterms com- 
pensating the imaginary parts arising from infrared renormalons and, in connection with 
the case of the ultraviolet renormalons, we check that a unique prescription can be given. 
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A short comment is in order about the use of the 1/Nf expansion. This kind of 
expansion is the natural framework for the study of renormalons in that it provides a 
graphic tool (see figures) for the investigation of these problems^. Moreover it allows a 
partial resummation of the perturbative series so that one can obtain informations on 
the whole series and calculate quantities, like the overall normalization constant in ( |1.4| ), 
which cannot be found at any finite order in perturbation theory [|12 |. 

The paper is organized as follows. In Sect. II we introduce the notations and study the 
problem of the first ultraviolet renormalon. We find the explicit form of the counterterms 
and we check that all singularities can be absorbed in such a way. Moreover we verify 
that the non local singularities which appear in the naive perturbation theory are indeed 
absent when considering the lagrangian with the counterterms. In Sect. Ill we perform 
similar calculations for the infrared renormalons; the main difference is that, in this case, 
the operators to be added to the lagrangian are multilocal. Sec. IV is devoted to some 
considerations about the meaning of the lagrangian with the imaginary counterterms. 
Finally, the Appendix contains some technical details about the calculations with the 
background field formalism. 



II. ULTRAVIOLET RENORMALONS 



Let us begin this section with the introduction of some definitions about the 1/Nf 
expansion in QED. 

Defining a rescaled coupling constant a = aNf, the 1/Nf expansion is generated, at 

— 1/2 

fixed a, by considering a factor of Nf for every interaction vertex and a factor of Nf 
for every closed fermion loop. Then the photon propagator is given at leading order by 
the sum of the bubble diagrams (fig.l). In the Landau gauge it reads: 

D, u {k) = - ^ K -f 9 ^ D{k) (2.1) 

where: 

D{k) = T+Tm*) < 2 ' 2 > 

and IIo(A:) is the single bubble diagram. In the following we are interested in the high 
momentum region of integration in the Feynman integrals; so we need the expression of 



After completion of this work, we came at vision of 25] in which the authors conclude that it 
is not possible in general to associate the ultraviolet renormalons to a well-defined set of graphs; 
this is indeed possible in the large Nf limit; on the other hand the calculations performed in 
our paper are at fixed order in 1/Nf. 
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Tlo(k) in the deep euclidean region. After renormalization in the MS scheme, we have (all 
fermions are assumed to have the same mass): 



n (fc) 



\ k 2 5 
log— - 



2 3 



2a m 2 

TT k 2 



+ 



(2.3) 



From ( |2.2| ) and the definition ( |1 . 1| ) one can easily derive the expression for the Borel 
transform, with respect to a, of the photon propagator; introducing an extra factor of a, 
we have: 

„2 \ k/Svr 



F(k,b) = (aD)(k,b) ~ e- 5b/97T [— 

fc>m \ \X 



2bm 2 „ ( 1 



(2.4) 



A. Imaginary part of Green functions 

At next-to- leading order in 1/Nf, the Green functions correspond to the graphs with 
one insertion of the chain of the bubbles (fig. 2-4) and present renormalon-type singu- 
larities. In fact the Feynman integrals, in which the photon propagator ( [2.2D appears 
in the integrand, have an ambiguity arising from the region of integration around the 
Landau pole (i.e. the pole in (|2.2|)) which occurr, for infrared free theories (in our case 



(3(a) = a 2 /3n + 0{1/Nf)) in the ultraviolet region: 

A2 = ^ e 5/3 e l//3 2 a (25) 

(we are always assuming the MS scheme). This ambiguity causes the Green functions to 
develop an imaginary part and is related to the leading factorial growth of the coefficients 
of the perturbative series in a. In fact it turns out that, after expanding (|2.2j ) in powers of 
a, the coefficient of order n receives a contribution proportional to n\ from the ultraviolet 
region of integration. These problems appear as poles on the real positive axis of the 
Borel transformed Green functions. 

In the following we intend to find the behaviour of the Borel transform of the vertex 
function, of the electron propagator and of the photon vacuum polarization near the first 
ultraviolet renormalon at b = 3tt. This will enable us to find the imaginary part of these 
functions and to study the connection with local gauge-invariant operators of dimension 
six. 

Let us begin with the Borel transform of the vertex function (fig. 2). In the euclidean 
space the term of order 1/Nf is given by: 

Ijt\ = 4tt r d 4 k -f p [p' — k — m^p — k — m]^ k p k a - k 2 5 pa 
e) N f J (2tt) 4 \(P' ~ k) 2 + m 2 ]\(p - k) 2 + m 2 } k 4 [,) 



where we have factorized out the electron charge e. Renormalization is taken into account 
by simply introducing the renormalized bubble. In fact, since we are working in the 
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Landau gauge, the vertex function needs no further renormalization and this results in 
the absence of a pole in b = in T M as written in ( |2.6|) . Moreover introducing the term of 
order l/k 2n in the expansion ( |2.4| ) of F(k, b) results in a pole in T M at 6 = 3ir(n + 1). 

After some manipulations with 7 matrices and performing the integral we obtain for 
the first pole: 

+Sp%p-4:(p^-p fM )(p' -p)} (2.7) 

From the Borel representation ( |1 . 1| ) we deduce that T M develops an imaginary part given 
by the residue at the pole in b = 3n: 

Imr M = 7reRes|^^e- b/a | (2.8) 

Then Q yields: 

Im T M = {[Aq 2 - 3p 2 - 3p' 2 - 9m 2 } lfl + 3p' lfl p - Aq^q] + O (J-^j (2.9) 

On the mass-shell this expression reduces to: 

Imr M = — -^g 2 7u (2.10) 
N f 2A 2 ^ ^ v ; 

hence there is no contribution of the first ultraviolet renormalon to the large order be- 
haviour of the perturbative series of the anomalous magnetic moment of the electron, in 
agreement with the asymptotic behaviour found in |J. 

From ( |2.9p we see that the imaginary parts arising from different renormalons are 
classified, at this order in 1/Nf, in terms of the renormalization group invariant (i.e. \x- 
independent) A-parameter which plays the role of a dimensional expansion parameter. It 
must be noted that this parameter encloses all the dependence of the calculation from the 
choice of the renormalization scheme. 

Similar considerations hold for the fermion propagator (fig. 3). Its Borel transform at 
order 1/Nf reads: 

w„ h) _ 47r f d ' k ^[(P + k)-m} lu k 2 5^-k^k v 

MP ' b) ~ W f J (2tt)4 (p + k) 2 + m 2 jfc* * ^ b) {2A1) 

Isolating the pole at b = 3n we have: 

1 / 1 \ b/37T e ~5b/9n (- g g 3 9 bm 3 ) 
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This yields an imaginary part given by: 



Im E(p) = ~ {- 9 -m*p - \p*p + ^m 3 } + O ) (2.13) 



It is straightforward to verify that (|2.9| ) and ( 2.13|) satisfy the Ward identity: 



g M Im ( -^) =ImS(p') -ImS(p) (2.14) 



e 



Finally, the imaginary part of the photon vacuum polarization can be deduced from 
Beneke ||24j| . Exploiting gauge-invariance for the proper photon prapagator we have: 



n F (g 2 ) = (^-9V) n (? 2 ) (2-15) 

with: 



n-Db + n^ + o^i (2.16) 



Il is the single bubble, while rix corresponds to the diagrams of fig.4. The result found 
by Beneke in the massless limit is: 



,2 \ b / 37r 

a J v ^ ' 6ir 2 \ fi 2 



Ul ^(Q 2 ,b) = --V 5 ^ ( % ) D(b/3n) (2.17) 



with: 



D(u) ~ \-r^— 2 + + 0(1) (2.18) 

3(1 — u)^ 9 1-w 

From this expression one can infer the imaginary part of the proper photon propagator: 

a q 2 ( q 2 5 \ 

Im ni^(g 2 ) = - - — {q^ ~ q 2 g^) ( - - J 

+ ^(g„ff„ - 9 2 V) + t) + O (p) (2.19) 

The non local singularity corresponding to the term proportional to logg 2 in the above 
formula could in principle destroy the connection with the gauge-invariant counterterms. 
Its role will be discussed in the following. 

B. Connection with local operators of dimension six 

Our basic result is the explicit value of the imaginary parts of the vertex function 
(eq.( [2.9| )), of the fermion propagator (eq. (|2.13| )) and of the photon vacuum polarization 



7 



(eg. ( p. 19] )) in correspondence of the first ultraviolet renormalon. We want to stress again 
that all our calculations have been performed at fixed order in 1/Nf and all the following 
considerations rely upon this kind of expansion in which renormalon-type problems have 
a more direct interpretation. 

Let us now turn to the main point of our investigation i.e. the connection with local 
gauge-invariant operators. On simply dimensional grounds and exploiting renormalization 
group arguments one can argue that the imaginary parts given by (f2.9|), ( 2.1^ ) and ( |2.19| ) 



are proportional to the insertion of operators of dimension six. We want to check if 
it is possible to give a prescription in terms of local gauge-invariant operators added as 
imaginary counterterms in the lagrangian. Let us study the structure of the counterterms. 

The term proportional to (<7 2 7 M — q^q) in ( |2.9|) can be seen as the insertion of the 
gauge-invariant operator (d^F ^ip^ip; hence we deduce the presence of the following 
counterterm in the lagrangian (the coefficient is uniqely determined): 

1 7fg _ 

AA = -i Wf ^(9,F^)^ (2.20) 

Moreover, from the term proportional to m 2 in ( |2.9| ) and ( |2.13| ) we see that the operators 
m 2 K~ 2 ^pdtl) and em 2 A~ 2 ijjAip combine in a gauge-invariant form and this corresponds to 
the following lagrangian counterterm: 

1 Qtt m 2 — 

A£ 2 = -i— ^^(id-eA)^ (2.21) 

The fundamental check involves the imaginary part of the photon vacuum polarization 
which contains a term proportional to logg 2 . In fact this term cannot be cancelled by 
the direct insertion of a local operator but ought to arise from the insertion of a local 
operator in a loop diagram. One can easily convince himself that the only operator that 
can do this job must have the form of ( 2.20Q . We then compute the insertion of AC\ on 



the vacuum polarization as given by the diagram of fig. 5. This diagram is of order 1/Nf 
since it contains a term 1/Nf from ACi, a term 1/Nf from the two vertices and a term 
Nf from the fermion loop. It gives the following contribution to the imaginary part of 
the vacuum polarization (in the MS scheme): 

1 a q 2 , 9 . / q 2 5\ . 

! ^3^-*V)(log^-3j (2.22) 

This contribution exactly cancels the logarithmic term in Imlli. 

To summarize, we have explicitly checked that, at next-to-leading order in 1/Nf, the 
imaginary parts due to the ultraviolet renormalon at b = 3tt can be cancelled by the 
contribution of local gauge-invariant operators of dimension six. Indeed it is possible to 
give a realization of this connection at lagrangian level. 
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III. INFRARED RENORMALONS 



In this section we deal with the so-called infrared renormalons i.e. singularities on 
the real axis of the Borel transform at bfa = —n which arise from the low momentum 
region of integration in the Feynman integrals. As already stated these singularities are 
on the negative real axis of the Borel variable for infrared free theories like QED or QCD 
in the limit of large number of flavours and give rise to a sign alternating Borel summable 
contribution to the perturbative series (see eq. (|1.4j )). 

On the contrary, for asymptotically free theories like QCD the infrared renormalons 
destroy the Borel summability of the perturbative series, a phenomenon usually related 
to the presence of non-perturbative effects, the condensates, which cause the physical 
vacuum to have a non trivial structure. Hence the study of infrared renormalons in QCD 
is of primary interest expecially for the comprehension of the operator product expansion 
which is usually introduced in order to parametrize non-perturbative effects [^J . 

Neverthless it is common wisdom that the structure of the infrared renormalons in 
QCD can be obtained from the study of QED, after the substitution of the value of fa 
of QCD. This means that the structure of the infrared renormalons is insensitive to the 
presence of non-perturbative contributions in the operator product expansion which in the 
case of QED is a pure reorganization of the perturbative series but depends only on the 
nature of the gauge-invariant operators of the theory, in agreement with the conjecture 



made by Parisi [20 



In the following we intend to perform an explicit check of the connection with gauge- 
invariant operators, in analogy with the calculations of the previous section concerning 
ultraviolet renormalons. The basic difference in the case of infrared renormalons is that 
this connection can be implemented by adding imaginary multilocal gauge-invariant op- 
erators to the lagrangian. These multilocal operators must be considered as the "dual" of 
the local operators appearing in the operator product expansion (for a rigorous definition 
201). Our aim is to find the form and the normalization of these operators and to 
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check the combinatoric. 

In analogy with the preceding section we deal with QED at next-to-leading order in 
1/iVj-. The only difference is that now we consider the massless case; then the proper 
photon propagator at leading order has the following form in the MS scheme: 



a / k 2 5\ 
n °« = -3?K?-3) <"> 



while the Borel transform of D(k) is: 

F(k, b) = e~ 5b/9n rjpj (3.2) 

The expression for the Borel transformed functions is the same as before. From ( |2.6|) and 



' L.2 \ b / 3n 
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[2J) we see that T M has infrared singularities for b = — 37rn; expanding around the first 
two poles we have: 



N f J [8 b + 37r p lfl p> 4 6 + 6tt 



2q 2 

PluV + 



+ 3 7 _ 5 v _ 5 M _v^_M. 

f p^ < f p^ ! ^' f p^' 



(3.3) 



This yields an imaginary part given by: 

Nf 8 p p' A*/ 4 p^/^ \p z p' z 

While for the ultraviolet renormalons the expansion parameter was A -2 , now it is A 2 ; 
in fact we are simulating ultraviolet free theories which present a Landau pole at small 
momenta (in ( |2.5| ) we must take (3 2 < 0) 

The imaginary part of the fermion propagator can be obtained in a similar way. We 
have: 

T N 1 9ttA 2 ^ 1 3ttA\ ^ /a , n . , 

ImE(rt = iv;T^-iv;T7 p+0(A) (3 - 5) 

Even (|3l|) and satisfy the Ward identity ( gig ). 

From p4j we deduce the following imaginary part for the photon vacuum polarization: 

Im n lMl/ = -j-^foq* ~ q 2 9,u) + 0(A 6 ) (3.6) 

The absence of the term of order A 2 reflects the absence of the first renormalon pole 
at b = 37r and this is a first evidence for the connection with gauge-invariant operators 

ESS 



We want to study how the imaginary parts ( |3.4| ), ( |3.5[ ) and (|3.6|) can be compensated 
by the contributions of multilocal operators. 

First of all it is evident that the terms of order A 2 in ( |3.4| ) and in ( |3.5| ) are both 
proportional to the insertion on the functions calculated at leading order (without the 
chain of bubbles) of the multilocal operator: 

ASl = ~ l W f li A2 / d * x ^W J di y^Mv) (3-7) 

It is pleasant that a unique choice of the normalization of the operator takes care of both 
the contributions. 
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Let us now turn to eg. (pTq) . When considering the operator product expansion for 
the vacuum polarization, the first non trivial operator to be considered is F^F^. It has 
dimension four and thus gives a contribution proportional to q~ 4 that should compensate 
the imaginary part of the vacuum polarization arising from the infrared renormalon at 



b^2 = —2 |22|| . We want to extract this contribution adding to the action the multilocal 
operator dual to F 2 : 

C 



AS 2 = t-A 4 J d 4 xF^{x) J d^F.M (3.1 



The constant C is fixed by the condition that the insertion of AS2 on Ho cancel the 
imaginary part of IIi. 

In general the insertion of AS2 on some correlation function can be computed with the 
Schwinger operator method in the so-called fixed-point gauge; the details of the calculation 
will be given in the appendix. The result for the insertion of AS2 on Hq^u is: 

Comparing this expression with eq.(|3.(j|), we find for the constant C the value: 



° = l < 3 - 10 > 

As a fundamental check for the consistence of the theory we try to extract the same 
constant C from the imaginary part of the other functions. In fact, while there is no 
contribution of AS2 on the fermion propagator (this is consistent with the form of ImE), 
one can compute the insertion of AS2 on at q = 0. Referring again to the appendix 
for the details of the calculation we find the following contribution: 



i e ^K^U-A P -4) (3.11) 



This must be confronted with the term proportional to A 4 at q = in ImT^: 



Hence the same coefficient C found in (|3.10|) eliminates this imaginary term. 

For the sake of completeness we have extracted the coefficient C directly from the 
imaginary part of (F^F^); from the diagram of fig. 6 we find for the Borel transform: 

(aF 2 ) ~ -f j-^-e-^l-A (3.13) 

which implies: 
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Im (F^F M ,) = -|-A 4 + 0(A 6 ) 

OCL 



(3.14) 



Moreover the insertion of AS2 on (F^F^) is simply given by C so that even this imaginary 
part is cancelled with the choice ( |3.10| ). 



IV. CONCLUSIONS 



We have shown in a non trivial context that renormalon-type singularities can be 
absorbed in an effective lagrangian of the form: 

C eff = C + iAC (4.1) 

AC is uniqely determined by the request of having real correlation functions in the eu- 
clidean region. Local counterterms correspond to ultraviolet renormalons while multilocal 
counterterms to infrared ones. 

The imaginary part in £ e // restores unitarity which is violated in the original theory 
for the presence of renormalons. 

The description in terms of an effective lagrangian makes sense only if the resulting 
theory doesn't present non local singularities. We have checked that this is the case for 
QED in the framework of the 1/Nf expansion. 

From the theoretical point of view it would be interesting to study the global structure 
of C e ff and to understand how the introduction of imaginary counterterms may influence 
the renormalization of the theory. 
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APPENDIX A: BACKGROUND FIELD TECHNIQUE 



We want to shortly review the background field formalism which enables to compute 



the correlation functions in an external field p6| , |27 |. 

The Schwinger method enables to enforce gauge-invariance at all stages of the calcula- 
tion; it is based on the introduction of a basis of eigenvectors of the coordinate operator: 



d 

{x\ p Ay) = ~ i ~faT S \ x -y) + eA^(x)5(x - y) 



(Al) 



In addition we perform our calculation in the so-called fixed-point gauge p8 

XpA^x) = 



(A2) 



In this gauge the vector potential is written in terms of the field strenght and of its 
covariant derivatives calculated in x = 0; in momentum space we have: 

A,(k) = J d 4 ze lkz A,(z) ~ -i^F pix (0)-^6^(k) + ... (A3) 

From this formulas it is easy to derive the expression for the fermion propagator in an 
external field: 

S(q) = I dx(x\-J—JO) = -- \eF pa q plal5 + ... (A4) 
J P+q IT 

where F Mi , is the dual of F pu . 

We want to find the contributions of (F 2 ) to the proper photon propagator and to the 
vertex function; for this purpose we must use the randomness of the field i.e.: 

(F^(0)F P(T (0)> = —{g w g va - g^g, P )(F a ^ aP ) (A5) 

The insertion of (F 2 ) on IT(g 2 ) has been computed many times in the literature (see for 
example [p3jP ; the expression for the photon propagator in an external field is given by 
(using (0) and flA5|)): 

IV(g) = ie 2 j d 4 xe^ x Ti{ lp S(x,0) lu S(0,x)} = (A6) 

The contribution of (F 2 ) to the vertex function at zero momentum transfer comes from 
the insertion of a single A p on every propagator (see fig. 7): 
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vp(p) 



-2e J 



d 4 ki d 4 k2 



A{k 2 ) 



P + K\ p + K\ P 



" v ^ ' (2tt)* (2tt)* £ + A* + A* 

From (|A3|) and (|A5|) we have, considering the one-particle irreducible function: 



Tf\p) = \e 



s<5> 
P 4 



which corresponds to eq.( |3.11| ) 
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FIGURE CAPTIONS 

Fig 1. The photon propagator at leading order. 

Fig 2. The vertex function at order 1/Nf. 

Fig 3. The fermion propagator at order 1/Nf. 

Fig 4. The photon vacuum polarization at order l/Nf. 

Fig 5. Insertion of ACi (the crossed circle) on the vacuum polarization. 

Fig 6. Graph contributing to the vacuum expectation value of F^ V F muv . 

Fig 7. Contribution of F 2 to the vertex function at zero momentum transfer. 
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